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ABSTRACT 

Like many young supernova remnants, SN 1006 exhibits what appear to be clumps of ejecta close to 
or protruding beyond the main blast wave. In this paper we examine 3 such protrusions along the east 
rim. They are semi-aligned with ejecta fingers behind the shock-front, and exhibit emission lines from 
O VII and O VIII. We first interpret them in the context of an upstream medium modified by the 
saturated nonresonant Bell instability which enhances the growth of Rayleigh- Taylor instabilities when 
advected postshock. We discuss their apparent periodicity if the spacing is determined by properties 
of the remnant or by a preferred size scale in the cosmic ray precursor. We also briefly discuss the 
alternative that these structures have an origin in the ejecta structure of the explosion itself. In 
this case the young evolutionary age of SN 1006 would imply density structure within the outermost 
layers of the explosion with potentially important implications for defiagration and detonation in 
thermonuclear supernova explosion models. 

Subject headings: acceleration of particles — cosmic rays — ISM: supernova remnants — shock waves 
— supernovae: individual(SN 1006) 



1. INTRODUCTION 

Supernova remnants (SNRs) have long been suspected 
as the sites of cosmic ray (CR) diffusive shock acceler- 
ation, where ^10% of the supernova energy would be 
required to be transmitted to CRs to satisfy observa- 
tions. Partial confirmatio n of this idea came with ASCA 
obser vations of SN 1006 (|Kovama et al.l 119951: iRevnoldsl 
|1996[) . which revealed NE and SW limbs dominated by 
X-ray synchrotron radiation from CR electrons. Direct 
evidence of CR ions, presumably the dominant CR com- 
ponent, has been harder to find, for the most part be- 
ing limited to the observation of 7-rays resulting from 
the rea ction p + p ^r" + ^ followed by tt" 2j 
decay. lEriksen et al.l ()2011[ ) interpreted newly discov- 
ered, regularly spaced stripes of non-thermal emission 
in Tycho's SNR as regions where the pre-shock medium 
has been significantly modified by ^ lO^^eV ions. Indi- 
rect evidence for CR ions comes from the separation be- 
tween the for ward shock and the contact discontinuity in 
youn g SNRs (IWarren et al.l [200l [C assam-Chena'i et al.l 
I2OO8I: iFerrand et al.l I20100 . For example, in SN 1006, 
the contact discontinuity is predicted to be less than 
0.86 of the forward shock radius, but is observed to 
be larger than that around almost all of the remnant 
(jCassam-Chenai et al.l 120081: iMiceh et al.] I2OO90 . This 
has been interpreted in SN 1006 and Tycho's SNR 
as bei ng due to energy losses to CRs at the forward 
shock (IWarren et all l2005t iCassam-Chenai et al.l I2007L 
[2008'; 'Mic eh et al.ll2009f ). which leads to a greater shock 
compression ratio and a thinner shell of shock ed CSM 
()Blondin fc Ellison! [2OOII: iFraschetti et aLllMot) . How- 
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ever, only in cases of extreme CR energy losses can a 
fraction of the ejecta in Rayleigh- Taylor (RT) instability 
approach or overtake the forward shock. 

Density inhom ogeneities in either the inter- 

stellar medium (lJun et al .' '1996 1) or the ejecta 
(|Wang fc Chevahed l200Tr 'Blondin eTaH l200l can 
also enhance the growth of RT instabilities potentially 
allowing them to reach the forward shock. Originally 
motivated to explain the existence of radi al magnetic 
fields in young SNRs pun fc Norman! [l996D . lJun et ahl 
(|1996l) explored the effect of turbulence induced as the 
SNR forward shock encounters density clumps (of a 
factor of 5) in the preshock medium. The vorticity 
induced by the interaction of the forward shock with 
the overdense clumps is advected postshock. Upon 
reaching the contact discontinuity, the vortices transfer 
their rotational energy to the RT fingers. This enhances 
the growth of the fingers such that ejecta clumps may 
protrude beyond the blast wave. In some specific cases 
it was found that twenty times larger density inhomo- 
geneities would be required in the eje c ta to achieve 
the same result (|Wang fc ChevalieH [20011: iBlondin et all 
I2OOI . 

Upstream density clumping may be intrinsic to the 
pre-shock medium, or induced by a CR shock precur- 
sor. In the scenario of nonresonant magnetic field ampli- 
fication (Bell 2004, 2005), density clumping by a factor 
of ~3 is produced upstream of the shock as the mag- 
netic field amplification saturates. Magnetic field am- 
plification, also taken by some as indirect evidence for 
proton acceleration at shocks, is suggested by the nar- 
rowness of r ims of synchrotron radiation at SNR shocks 
(iVink fc La ming 2003; Long et al. 2003; Yamazaki et ^ 



2004!: IWarren et al.. .2005; .Cassam-Chenai et ah ,2007 



Volk et all 12001 ? Modeling the widths allows one to 
estimate the magnetic field to be in the range 100 - 500 
fiG^ far larger than can be attained by the typical com- 
pression of ambient magnetic field. Fields of this magni- 
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tude may speed up the CR acceleration process, allowing 
SNRs to generate CRs up to and beyond the "knee" in 
the CR spectrum, at about lO^^ - lO^^ eV. 

Of all known SNRs, SN 1006 offers the most promis- 
ing test of the action of CRs to produce density clump- 
ing ahead of the forward shock. It is the likely rem- 
nant of a Type la supernova, located 500 pc above the 
galactic plane and expanding into low density interstel- 
lar mediu m. A H.E.S.S. detection of SN 1006 has been 
reported (jNaumann-Godo et al.l l2008| ) , consistent with 
an ambient gas density of 0.05 cm~'^ (jAcero et al. 2007). 
Proper motions hav e been measured in the NW in the 
optical, 0.28" yr"! (|Winkler et al.ll2003f ) corresponding, 
at a dis tance of 2.2 kpc, with a shock velocity of 2900 
km s~^ (|Ghavamian et al.l[2002l ). I n the NE, 0.48" yr~^ 
expansion was measured in X-rays (jKatsuda et al.|[2009l) 
giving 5000 km s~^, indicating higher density in the NW 
than elsewhere. Hamilton et al. (1997; 2007) derive a re- 
verse shock velocity of 2700 km s^^, and determine the 
expansion velocity of ejecta entering the reverse shock to 
be 7000 km s~^, from HST observations. 

In this paper, we examine Chandra images and discuss 
structures seen to protrude beyond the eastern portion of 
the shock front. Ejecta clumps with similar periodicities 
are seen behind or near the blast wave elsewhere in SN 
1006, most prominently around the southern rim. We 
chose to study these three features because they extend 
past the main shock and have deep enough observations 
for minimal spectral analysis. A more definitive study 
of all periodic structures in SN1006 is defered to future 
work, thus numerical results reported herein may not be 
apphcable for all of SN1006. 




Fig. 1.— From left to right, the 2000 and 2008 Chandra obser- 
vations of the eastern portion of SN1006 from 0.2 to 5.0 keV and 
a merged image whe re t he 2008 ima ge wa s shifted by the proper 
motion measured bv lKatsuda et^aP ||2009I) . All three panels were 
binned by 4 pixels, smoothed and use square root scaling to show 
both the faint and brighter features. The overlaid regions show 
where spectra were extracted from the 2008 observation. 



2. OBSERVATIONS 

The eastern portion of SN1006 has been imaged mul- 
tiple times with Chandra and XMM, and the features 
under discussion here are evident in all observations. 
However for the sake of higher statistics we focus only 
on the two deepest Chandra observations in 2000 July 
10, (Obsid 732) and 2008 June 24 (Obsid 9107). As 



TABLE 1 

Size scales of upstream features 





angular 


linear 




(arcseconds) 


(IQi* cm)" 


azimuthal length 


83, 41, 46'' 


2.7, 1.5, 1.4 


radial extent 


25, 20, 10 


0.8, 0.7, 0.3 


spacing*^ 


145, 125, 125 


4.8, 4.1, 4.1 



^ Assuming 2.2 kpc distance to SN 1006 
^ Three features, from south to north 

Measured from "leading edge" to leading edge of each feature as 
opposed to the size of the darker regions between features. The 
northernmost size is measured from the feature marked with a 
white arrow in Figure 1. 

in IKatsuda et al.l (|2009D we exclude times of high back- 
ground from the 2000 observation but found no flares in 
the 2008 observation. To further reduce the background 
we consider only the band in which the emission from 
SN1006 dominates, from 0.2 to 5.0 keV for imaging, and 
0.3 to 6.0 kcV for spectral analysis. 

Figure 1 shows the region of interest in each long ob- 
servation and a merged image where the 2008 observa- 
tion was shifted 3'.'85 along a line 7.5° S of W to match 
the proper motion found in IKatsuda et al.l ()2009D . The 
features outside of the main shock are highlighted with 
green polygons that mark the spectral extraction for the 
2008 observation. (The 2000 spectra were extracted from 
proper- motion-shifted regions). The southernmost fea- 
ture extends below the spectral extraction region, inter- 
secting the bright shock below the chip gap. The co- 
ordinates of the polygons are the same in both frames, 
such that one can clearly see the movement of the outer 
limits of the features from 2000 to 2008. The proper 
motion of the features appears comparable to the main 
bright shock, but is difficult to determine exactly given 
their low surface brightness. A fourth protrusion to the 
north that was too narrow to extract a spectrum from is 
marked in the third panel with a white arrow. Similar 
features can be seen further north. They do not extend 
so far in advance of the forward shock, and have no ob- 
vious connection with structures extending post shock. 

In Table 1 we present the size scales of the upstream 
features. The spacing between the three features corre- 
sponds to a wavelength of ~ 4 x 10^* cm at the distance 
of SN 1006 and assuming all are in the plane of the sky. 
Looking to the remnant interior one can see clumpy ex- 
tended fingers of brighter emission with a similar spacing 
underlined with dashes in Figure 1. The most obvious 
of these extends from behind the southernmost feature, 
but others are suggested behind each feature. The peak 
brightness of the protrusions is between 15% and 45% of 
the peak brightness of the main shock front immediately 
behind it. 

Figure 2 shows the combined spectrum of all three fea- 
tures ahead of the shock in the 2000 data. A background 
region from the same chip that excluded the features was 
chosen and modeled rather than subtracted and the spec- 
tra were binned to allow the use of statistics. The 
total counts in the 2000 data between 0.3 and 6.0 keV 
were 1650 and 1033 for the 3 features and the background 
region respectively. Table 2 summarizes the best fit val- 
ues and 3cr confidence intervals for a simultaneous fit of 
the 2000 and 2008 data. An absorbed power-law model 
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Fig. 2. — The combined spectrum of all three regions "ahead" 
of the shock from the 2000 observation, fitted with an absorbed 
power-law model plus two Gaussian lines for the H-like and He- 
like O. 



TABLE 2 

Spectral model fits for the J2000 and J2008 data 
Parameter Powerlaw + Lines Miceli SE, Powerlaw -t- NEI 



Nh (cm-2) 


5.8+« X 10^0 


7 X 10^" 


7 


r, 7+0.6 
^■'-0.4 


r, 7 + 0.3 
^■'-0.5 


Line (eV) 


56111? 




Line (eV) 






kT (keV) 




0.57+0-6 


net (cm~^ s) 




5.7tl\ X 108 


(dof) 


200.5(197) 


210 (199) 



with two additional Gaussian lines fits the data well, with 
an appropriate Njj for SN 1006, and a power-law in- 
dex consistent with fits to a nearby region of the bright 
shock. The addition of two Gaussian lines for He-like 
and H-like O is supported by an Ftest at the 99% sig- 
nificance leve l. Alt ernatively, following the SE model of 
iMiceli etTaLl (|2009f ). the data are fit equally well with 
a synchrotron and non-equilibrium ionization thermal 
component. In this fit we allowed only the normaliza- 
tions, powerlaw index 7, temperature kT, and ionization 
timescale net to vary. Their analysis yielded a super-solar 
oxygen abundance which matches our spectrum well, the 
only difference being the best-fit ionization timescale for 
the 3 features is nominally lower than that of the shocked 
region as a whole (at la, but not at 3cr). 

3. DISCUSSION 

3.1. Overview 

The most relevant observed characteristics of the three 
pre-shock features are the generic size-scale of ~ 4 x 
lO^^cm, a continuum that is consistent with the syn- 
chrotron emission from the bright shock, and the pres- 
ence of He-like and H-like O lines that imply an ion- 
ization age of ~ 5.7t-1\ x 10^ cm^^ s. This ionization 
age, comparable to the remnant lifetime given the ambi- 
ent density of SN 1006, argues against an interpretation 
of the features as the denser portions of a CR precur- 
sor. The likely supersolar abundances with a synchrotron 
dominated continuum, relative brightness of the preshock 
features, and small radius of curvature also support the 
idea that these are ejecta fingers that have breached the 



shock front. 

Hereafter, we discuss possible origins for the spacing 
of the pre-shock features. (1) Anomalous viscosity deter- 
mines the preferred RT wavelength. (2) The saturated 
state of nonresonant magnetic field amplification creates 
cavities and higher density regions at a particular scale. 
(3) The ejecta fingers originate in the explosion itself. 

3.2. Ejecta Fingers Spaced by Anomalous Viscosity? 

If the preshock features seen are the extension of the 
ejecta fingers themselves (jCassam- Chenaj- et al.l [20M) 
then the spacing between them may reflect the wave- 
length of the R-T instability. In the absence of a peri- 
odic perturbation, this depends on the kinematic viscos- 
ity (the product of the sound speed and largest eddy size: 

^ — ^sound^eddy) • 

In stationary incompressible media the Rayleigh- 
Taylor linear growth rate is approximately T = 
•y/ v'^k'^ + gk — vk^ which has a maximum at mode 

wavenumber k = (g/i^-^) ^^'^ /2 (iPlesset fc Whippli 

\1974} . where g ~ 0.0049 cm s'^ for SN 1006 is the decel- 
eration o f the plasma (estimated f rom the SNR dynam- 
ics using l^uelove fc McKeel[l999D . and the ejecta den- 
sity is assumed much greater than that of the shocked 
ISM. Exponent ial growth goes over to p ower law growth 
as displayed in lBlondin fc EllisonI ()2001[ ) when expansion 
is included. Taking the observed wavelength of 4 x 10^^ 
cm, we find v ^ 10^^ cm^s~^ and a maximum linear 
growth rate g'^/^/2v^/^ ~ 6 x 10"" s'^. This is slow, 
allowing only 2 e- folding times during the 1000 yr age 
of SN 1006. Additional vorticity is probably needed to 
speed up the process. 

We can compare this to an estimate of the viscosity 
from the product of the sound speed in the shocked inter- 
stellar medium (about 4 x 10^ cm s~^) and a length scale 
determined by the separation between the forward shock 
and contact discontinuity (about 1 pc) to find ^ 10^^ 
cm^s~^. Thus either the viscosity is anomalous by a fac- 
tor of ^ lO^'^, or the wavelength of the instability is 
determined by different physics, i.e. the CRs. 

3.3. Clump Spacing Determined by Cosmic Rays? 

We explore the case that the ejecta finger spacing is 
determined by structure in the CR precursor where non- 
resona nt magnetic field amplification has reached satu- 
ration (lBeil[200l l2005l) . Structure arises because the 
highest energy CRs, which remain unmagnetized inhabit 
cavities in the upstream medium which are evacuated by 
the amplified magnetic field. The plasma and magne- 
tized CRs are swept into cylindrical "walls" surrounding 
the cavities where the originally quasi-parallel shock has 
now become quasi-perpendicular. Interpreting the ob- 
served structures as ejecta fingers protruding into the 
cavities, the cavity radius, observed here to be around 
55", or 2 X 10^* cm at the 2.2 kpc distance of SN 1006 is 
similar to the unmagnetized CR gyroradius r^, so 



7mc 



2 X lO^^cm ^ B - 2 



(1) 



where 7 is the Lorentz factor for a high energy unmagne- 
tized CR and B is the amplified magnetic field. The CR 
rest mass and charge are m and q, and c is the speed of 
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light. To estimate B and 7, we consider a c ondition fo r 
the existence of the nonresonant instabihty ()BeUll200"^ 



< h < 



JcrB ATTUcFi^qvs 



pcv\ 



B < ATTmVsCncR 
(2) 

where = B^/Anp and Jcr = ncRqVg/l is the CR cur- 
rent in terms of the total number density of CRs ncR, 
(of which the unmagnetized fraction nj^^ = ucr/j), the 
charge q, and the shock velocity Vs at which CRs are 
assumed to stream^ For SN 1006 with Vs — 5000 km 
s-\ < 3 X IQ-^ncR. The total CR number den- 
sity in terms of 77 = Pcr/pv^, the ratio of CR pres- 
sure to shock ram pressure, with Pcr. — Ecr/S = 
^pZT fcRpdirp^dp/i, is 



ncR 



In 7C^ 



(3) 



where rii ~ 0.05 cm^'^ is the ion density in the upstream 
medium. Hence ^ lO"^/;, or with 77 = 0.1, B ^ lO^G 
then 7 ~ 6 X 10^. At a position further to the north, 
iMorlino et"all (|2010f ) infer a value of 77 ~ 0.29, which 
results in higher values oi B ^ ITuG and 7 ^ 10^. Both 
these values imply E = 10^^ — 10 eV, which is past the 
"knee" in the CR spectrum. 

3.4. Ejecta Fingers Associated with SN Explosion? 

The possibility remains that the ejecta structure inter- 
preted as due to RT instability above is in fact related 
to the Type la explo sion, as in e.g. the 0-DDT i nodel o f 
iMaeda et all (|2010a') illustrated in lMaeda et"all (j2010b[ ). 
In other SN la remnants, such as Tycho, similar clumps 
to the ones considered here show the clear presence of Fe 
ejecta and as a result have bee n interpreted in this way. 
Delayed detonation models (e.g lGamezo et al]|2004[) ap- 
pear t o be the most plausib le explosion model for a Type 
la SN (jBadenes et al.|[2006l see e.g.) . Ashes burnt by the 
deflagration m ay exhibit structure (|Maeda et al.ll2010al : 
[Jordan et al.l[2008. ). whereas large density gradients ap- 
pear to_be_ateentfirandert wave 
(e.g IMaeda et al.ll2010at iMeakin et all I2009D . Unburnt 
ejecta exterior to these regions may survive due to ex- 
pansion ahead of the deflagration front, with the den- 
sity diminishing sufficiently to allow them also to escape 
burning in any subsequent detonation as well. 

In SN 1006, both spectrally and dynamically, it ap- 
pears that the reverse shock has only recently encoun- 
tered regions that were subject to thermonuclear burn- 
ing. Suzaku recently detected Fe K emission with a low 
ionization age consistent with the Fe ejecta being more 
recently shocked than the other ejecta (Yamaguchi et al. 
2008). The ejecta expansion velocity at the reverse shock 
in the observations of Hamilton et al. (1997; 2007) is 
close to the predicted outer extent of Fe produced in the 
C-DDT and C-DEF models of Maeda et al. (2010a), and 
also that in the various models tabulated by Badenes et 

® We are taking the CR distribution function fcR = 
nCRPinj/^T^p'^, ncR = fp^'^'° fcR'i-n^P^dp with p^ax » Pinj, 

and n'p^ = /p^"."" fcR^'^p'^dp ~ ncRPi„j/Pmin ^ ncRh- 
The injection momentum is Pinj ~ mc, so that 7i„j ^ 1, and 
Pmin — 'ymc^ is the minimum momentum of the unmagnetized 
CRs. 



al. (2005). Hamilton et al. (1997) comment that the ISM 
density at the back of the SNR is likely to be lower than 
elsewhere, and so the reverse shock in other regions of 
the remnant may have progressed further into the ejecta, 
shocking some of the Fe as is now seen by Suzaku (Ya- 
maguchi et al. 2008). By comparison, in Tycho's SNR, 
Badenes et al. (2006) find that the reverse shock has 
reached all but the last 0.4 Msun of ejecta. Thus these 
la SNRs are at significantly different evolutionary stages 
in terms of the location of the reverse shock in the ejecta. 
An ejecta origin for the protrusions in SN 1006 would 
have different implications for the locus of deflagration 
and detonation burning than for similar clumps in Ty- 
cho, in terms of the location of deflagration ashes and 
whether their structure can survive a detonation wave, 
as in the asymmetrical 0-DDT model of Maeda et al. 
(2010a). 

4. CONCLUSIONS 

The structures observed exterior to the blast wave 
along the eastern limb of SN 1006 are likely to be metal- 
rich ejecta based on their probable supersolar abun- 
dances of O. We beli e ve th at a mechanism combining 
the ideas of lJun et all 11996') and lBelll (|2004[ ) is one way 
to bring the ejecta so close to the blast wave. Upstream 
density inhomogeneities of only a factor of 5 in the ISM 
could allow some RT fingers to overtake the blast wave, 
with the cavities and density enhancements that are a 
natural consequence of CR magnetic field amplification. 

We contend that prior explanations for ejecta over- 
taking the shock are insufficient if acting alone. In- 
voking enhanced s h ock co mpression following CR losses, 
iBlondin fc EUisonl H2001f ), needed post-shock compres- 
sions by a factor of 21 to obtain RT fingers that actually 
overtook the forward shock. This would imply implau- 
sibly high energy losses (77 % of shock ram p ressure lost 
as CRs, from equation 2 of iLiang fc Keiltvl I2000D . The 
recent more conservative simulations of iFraschetti et al.l 
(|2OT0f) do not show RT fin gers that exceed the blast 
wave radius, but do show that 2% of the ejecta mate- 
rial reaches with in 95% of the forward shock radius (see 
also lmngllMl . 

If the protrusions are RT fingers, we can interpret the 
size scale of the pre-shock features in terms of properties 
of the SNR in this sector of the shell. If the size scale is 
determined by the properties of the SNR plasma, assum- 
ing that the CR precursor itself establishes no preferred 
wavenumber for RT growth, then the viscosity would 
need to be anomalous by a factor of ~ 10~^. If ins tead 
the CR pr ecursor develops a preferred size scale (|Belll 
12004 I2005f) . and this determines the wavelength of RT 
growth, the spacing may be interpreted in terms of CR 
parameters to give B ^ 10/iG and 7 ~ 6 x 10®. This 
implies CR acceleration up to the knee in SN 1006. 

The other tantalizing possibility is that these struc- 
tures have an origin in the ejecta. If deeper observations 
revealed Fe as is seen in some ejecta clumps in Tycho's 
SNR, the dumpiness, Fe, and low evolutionary age of SN 
1006 would imply that defiagration ashes survived near 
the exterior of the exploding white dwarf. 
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